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Abstract 

We present a scanning tunneling microscopy (STM) study of native defects in graphene islands grown 
by ultra-high vacuum (UHV) decomposition of ethylene on Cu(l 11). We characterize these defects through 
a survey of their apparent heights, atomic-resolution imaging, and detailed tunneling spectroscopy. Bright 
defects that occur only in graphene regions are identified as C site point defects in the graphene lattice and 
are most likely single C vacancies. Dark defect types are observed in both graphene and Cu regions, and 
are likely point defects in the Cu surface. We also present data showing the importance of bias and tip 
termination to the appearance of the defects in STM images and the ability to achieve atomic resolution. 
Finally, we present tunneling spectroscopy measurements probing the influence of point defects on the local 
electronic landscape of graphene islands. 
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I. INTRODUCTION 


For logic applications looking to take advantage of graphene’s high mobility eharge carriers, 
much recent work has focused on overeoming two main challenges: incorporation of intrinsic 
dopants and introduction of a band gap in the electronic structure. Both challenges ean be ap¬ 
proached by manipulating graphene’s crystal structure. Compared to bulk, two-dimensional mate¬ 
rials offer a unique advantage since the eleetronic structure of the entire material can be modified 
through functionalization of its surface. This functionalization can proceed from adatom binding, 
the ereation of lattiee defects, or the substitution of atoms at C sitesHydrogen functionalization 
of graphene (to create graphane) shows evidenee of opening a band gap,® as was predieted by 
density functional theory calculations.® Sparse single vaeancies are shown to be effective sites for 
introducing dopants while maintaining reasonable transport properties.Di-vaeancies are also 
considered attractive sites for substitutional doping due to the size and band structure effeets of 
eertain potential dopants (sueh as Al, P, As, Ga, In, Pt, and Co), but experimental work is still 
limitedMost experimental work so far has involved boron and nitrogen substitutions at sin¬ 
gle C sites.nSQ^ These are either incorporated into the lattiee during ehemical vapor deposition 
growthpmi or afterward by ion irradiation and annealing and provide p- and n-type doping, re¬ 
spectively. A common challenge in functionalizing graphene through defeet engineering is in the 
defects’ relatively high formation energy and their tendeney to oluster,®i^whioh further motivates 
study in defeet ereation and eontrol. 

Defects in graphene are also of fundamental interest. Single vaeancies were predicted to have 
magnetie moments,!^ and there is now strong evidenee for such.ISJUl Point defect engineering 
may also enable control of plasmons on sub-wavelength seales.E^From a structural point of view, 
graphene and other 2D erystals provide new opportunities to directly correlate defeets and macro- 
seopie properties sueh as transport. For example, Qi and eoworkerd^SI used TEM and transport 
measurements to directly correlate defect healing at the edges of graphene nanoribbons with in¬ 
creased conduetivity. Prior TEM studies imaged the formation and mobility of graphene defeets in 
real-time,l^ studied the insertion of interstitials,!^ and probed doping via vacaneies ereated with ion 
bombardment.® Because of their relatively high formation energy,l22l graphene and graphite defect 
studies are typically limited to defects that are created by ion bombardment and irradiation. ®^32lll2l 
As a stochastie process, the defect density and elustering produced by ion bombardment provide 
limited eontrol. Native defects created during growth may allow for better control, but graphene 
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grown in near-equilibrium eonditions (sueh as CVD) is unlikely to exhibit native lattiee defeets, 
exeept at grain boundaries where independently nueleated growth regions meet. In one study, 
defeetive graphene was aehieved at intermediate phases during growth in UHV; the vaeaney eon- 
eentration eould then be eontrolled by healing of vaeaneies during subsequent thermal annealing.l^ 
Here we report a seanning tunneling mieroseopy study of native graphene defeets that form 
during ultrahigh vaeuum graphene growth on Cu(l 11). The growth proeeeds by ethylene deeom- 
position on a hot Cu(lll) erystal in short (5-10 minute) eyeles.EHSJxhe UHV-grown graphene is 
then studied without exposure to air by in situ transfer into the STM ehamber. We eondueted a sur¬ 
vey of over 100 point defeets on 15 different graphene regions from 3 similar graphene growths. 
Atomie resolution imaging helps us identify bright defeets as graphene lattiee point defeets, likely 
single vaeaneies. We present evidenee showing that dark defeets are most likely point defeets 
in the underlying Cu lattiee. The imaging eontrast and resolution are diseussed with respeet to 
bias voltage and the STM tip’s atomie-seale termination. We also present dl/dV speetroseopy 
probing the loeal density of states (LDOS) over point defeets, whieh reveals a relatively limited 
eontribution to the loeal eleetronie landseape. 


II. METHODS 

Reprodueible, epitaxial growth of graphene islands on Cu(lll) is aehieved by thermal de- 
eomposition of ethylene in ultra-high vaeuum (UHV) eonditions. This method produees pristine 
graphene that ean be studied by STM without exposure to air. We first prepare a elean Cu(lll) 
substrate by repeated eyeles of Ar+ sputtering and annealing at 600°C in 5-minute intervals to 
remove surfaee eontamination. Graphene is then grown by introdueing ethylene gas at a pressure 
of 2 X 10^® mbar into the preparation ehamber while eyeling the sample temperature from room 
temperature to 950°C 2-4 times. Eaeh temperature eyele is performed in 5 minutes and we pump 
the ethylene out of the ehamber when the sample temperature is at ~500°C during the eooldown 
of the last eyele. This proeedure leverages the enhaneed stieking of ethylene on Cu at room tem¬ 
perature, while promoting the the thermal deeomposition and self-assembly of graphene at high 
temperature. The sample is immediately eooled to 100 K and transferred into the STM ehamber 
where it further eools to 5 K in ~6 hours. 

The resulting graphene islands on Cu(l 11) are studied with a CreaTee STM at 5 K. Using STM, 
we ean probe topographie and speetroseopie information at an atomie seale. Images presented here 
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were obtained by measuring small changes in height of the STM tip to maintain constant current 
while scanning across the sample surface. A cut Ptir tip was used for all of these measurements. 
For spectroscopy measurements, the tunneling current was recorded while sweeping the bias volt¬ 
age at a tip height fixed by the initial tunneling conditions. By adding a small modulation of 5-10 
mVpp at --1 kHz to the bias voltage, a lock-in amplifier simultaneously measures the differential 
conductance, dFdV, proportional to the convolved tip and sample local density of states (LDOS). 
Spatial maps of dFdV at a specific voltage can be collected simultaneously with topographic im¬ 
ages by recording the lock-in signal pixel-by-pixel during scanning. 


III. RESULTS 

Five representative examples of graphene islands exhibiting native defects are displayed in 
Figure 1. Graphene is most often observed near step edges, and can grow continuously over 
terraces, as shown by the continuous Moire pattern over three terraces of Fig. lb. Islands often 
grow with a regular, hexagonal shape as shown in Fig. Ic. A range of island sizes was observed, 
with an average area of 940 nm^. The point defects in the graphene regions can be classified 
as either bright or dark by their contrast in STM images (c.f., yellow and blue arrows in Fig. 
la respectively). We surveyed 93 bright and 29 dark defects on 15 graphene islands from three 
similar, but independent growth runs. Additional features are sometimes seen, such as dark clusters 
of defects (Fig lb, red arrow) and irregular lines (Fig le, green arrow), which are found to be grain 
boundaries between graphene regions when imaged with atomic resolution. 

On graphene, we find an average of one dark defect per 450 nm^. A much higher density of dark 
point defects is seen on the Cu (e.g.. Fig. lb), and are commonly attributed to surface contaminants 
such as CCPSor O.ISl Auger spectroscopy of the sample after growth (not shown) indicates traces 
of S and O which may segregate from the bulk to the surface under the growth conditions. We also 
expect an amount of unreacted carbon from the ethylene decomposition. The dark point defects 
observed in the graphene regions have similar apparent depths and lateral extent as those observed 
on the neighboring Cu(l 11) substrate. They are also consistent with those observed for graphene 
on Cu(lll) by Martinez-Galera et al.,1^ which were attributed to defects in the underlying Cu 
surface. We will revisit this assignment below when we discuss atomic resolution images. 

The more common bright defects are only found on the graphene regions with a density of 
approximately one defect per 150 nm^, and will be the primary focus of this report. STM does not 
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have elemental speeifieity a priori, so indireet measurements are needed to gain further insight into 
the origins of the defeet. For example, we earefully eompared the apparent heights of the defeets in 
STM images. The apparent heights of the bright defeets vary from 8 pm to 21 pm under identieal 
tunneling eonditions (1 V, 0.2 nA), though they do not separate into elearly distinet groups. We 
attribute this range of measured apparent height to slight differenees in tip termination, whieh 
beeome more important for small apparent heights like those measured here. To eontrol for tip- 
dependent variations, we eompare apparent heights of defeets within the same image and imaged 
with the same tip termination. Figure 2a shows an example of a small hexagonal graphene island 
with 11 bright defeets and 3 dark defeets. This island is one of the most defeetive graphene regions 
we observed. Representative line profiles over a bright and dark defeet are shown in Figs. 2b and 
2e, respeetively. The apparent heights of the 11 bright defeets span a range with an average and 
standard deviation of 17 ± 3 pm. This is a mueh smaller range than in our larger survey, and 
suggests that there is only a single origin for the bright defeets we observe. 

To gain further insight into the atomie origins of the defeets, we eompared STM images under 
varying bias voltages (e.g., IV, 0.05 V, and -1 V in Figures 2d-f). The bright defeets in Figures 
2d and 2f are similarly imaged as eireularly symmetrie protrusions, though the apparent height 
ehanges from 21 pm to 26 pm respeetively. At low bias (Fig. 2e), the bright defeets only have an 
apparent height of ~13 pm, and are largely obseured by seattering of Cu(lll) surfaee state elee- 
trons. These eleetrons eomprise a nearly free 2D eleetron gas (2DEG) at the surfaee.l^ Reeently, 
seattering of these surfaee state eleetrons was used to infer the eharge state of adatoms on ultra- 
thin films of NaCl grown on Cu(lll).l2SlWe observe no additional seattering of the surfaee state 
eleetrons due to the bright defeets, whieh suggests that these defeets are eonfined to the graphene 
layer, and have little or no effeetive eharge. In eontrast, the dark defeets remain elearly visible in 
STM images at low bias, and we do observe eireularly-symmetrie seattering patterns of the surfaee 
state eleetrons. 

To probe the graphene lattiee in the vieinity of the defeets, we first note that atomie resolu¬ 
tion depends strongly on the termination speeies of the STM tip, whieh is not known under most 
eireumstanees. Although the tip is Ptir in the bulk, Cu or other adsorbates ean terminate the tip 
due to the repeated eontaet with the surfaee neeessary to optimize imaging. The tip termination 
ean be ehanged intentionally by eontaeting the sample surfaee and pieking up surfaee atom(s), or 
by deliberately approaehing e.g., adsorbed moleeules sueh as CO under appropriate eonditions.l^ 
Moleeule-terminated tips often provide sharper resolution than metal-terminated tips due to intro- 
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duction of additional chemical contrast. One of our motivations for using only sub-monolayer cov¬ 
erages of graphene, is that we have nearby clean Cu(l 11) that can be examined at the same time. 
In particular, we can look for anomalous features in tunneling spectroscopy such as prominent 
peaks that are often attributed to e.g., molecular orbitals. Furthermore, our experience suggests 
that atomic resolution on the close-packed Cu(l 11) surface is not achieved with metal-terminated 
tips under typical tunneling conditions; thus the appearance of the Cu lattice in STM imaging is a 
proxy for some kind of molecule-terminated tip. Because the tip termination often changes during 
the course of the experiments, we carefully look for changes in position or contrast during repeated 
imaging. 

STM images of an area imaged under the same tunneling conditions, but with three different tip 
terminations are shown in Fig. 3a-c. Tunneling spectroscopy and images on Cu suggest that Tip A 
was terminated with a metal atom (probably Cu), Tip B with an unknown molecule, and Tip C with 
a CO molecule. With Tip A, which is most common, the two defects in the image appear as bright, 
circularly symmetric protrusions. The same defects appear with atomic resolution when imaged 
with the molecule-terminated Tips B and C. For both tip terminations, the honeycomb lattice of 
graphene is clearly resolved and the defects are imaged as protrusions with trigonal symmetry. 

Additional substructure is shown in Figures 3d-e, which are STM images of the same defects 
with Tip B, but under different tunneling conditions (lower current). Since the tunneling con¬ 
ditions control both the tip-surface separation (via the current feedback loop) and the energy of 
tunneling electrons (via the applied bias), the resolution resulting from a molecule-terminated tip 
can depend strongly on these parameters. With this tip, we were able to tune the tunneling param¬ 
eters to achieve highly-resolved images showing the surrounding graphene lattice and the internal 
structure of the bright defect. The graphene honeycomb lattice is clearly resolved around the de¬ 
fects, where the centers of the carbon rings image as protrusions. By overlaying the graphene 
lattice on the images, we find that the defect is centered on a single carbon site, and has perfect 
trigonal symmetry extending from it. Fig. 3e shows that the two defects, which are oriented in the 
same direction with respect to the graphene lattice, occur on the same sublattice site (blue dots). 
Different orientations of the defects are also observed. Figure 3f shows an atomic resolution image 
of a different graphene region imaged with a different tip (Tip D), that appears similar to Tip C. 
Three triangular defects appear in Fig. 3f. Two form upward-pointing triangles and the third forms 
a downward-pointing triangle. Here, the honeycomb centers image as depressions, and we again 
overlay the graphene lattice to determine that the two defect orientations correspond to opposite 
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sublattice sites (lower insets). 

Overall, the defeet symmetry, orientation, and sublattiee assignments agree with theoretieal 
expeetations for a single vaeaney.l^SESI The three-fold seattering pattern from the defeets arises be- 
eause of suppressed intervalley seattering in graphene and refleets the symmetry of the sublattiee. 
Thus one would expeet that defeets on different sublattiees would be rotated by 60°, as observed. 
We eompare our images of single defeets to simulated STM images of defeets ealeulated using 
density funetional theory (DFT)PS1 and a eombined tight-binding approximation/DFT methoc^ 
for a freestanding graphene sheet. Perfeet trigonal symmetry is expeeted for an unreeonstrueted 
vaeaney.ES! However, a reeonstrueted form of the vaeaney, where the C-C interatomie distanee 
aeross the vaeaney shortens to reduee its overall energy, is expeeted. Simulated STM images of 
the reeonstrueted vaeaney show bright triangular features that laek three-fold symmetry,l22142l un¬ 
like what we observe. This suggests that the underlying Cu substrate stabilizes the symmetrie, 
unreeonstrueted form of the defeet. The alignment of the triangular features with respeet to the 
graphene lattiee imaged at low voltage is eonsistent with simulated images of the reeonstrueted 
vaeaney reported by Yazyev and eoworkers.l^ However, we note that the region of Fig. 3e imaged 
with the same tip but at IV (not shown) shows instead an orientation rotated by 60°, eonsistent 
with the simulations of Amara and eoworkers.l^SlThis differenee highlights the importanee of the 
bias voltage in interpreting experimental and simulated STM images. 

We note that nitrogen is a eommon residual gas in UHV, and one eould assign the defeets 
we observe to substitutional nitrogen instead. In faet, STM images of nitrogerPI substitutions 
also show the expeeted trigonal symmetry.EDSl However, the experiments and DFT simulations of 
these substitutions show image eontrast with symmetrie bright spots at the three nearest-neighbor 
sites to the N atom.E^ This differenee and their larger apparent heights (~80 prr|i^ makes them 
elearly distinguishable from single vaeaneies and our observations. In our system, the most likely 
alternative to a single vaeaney is a hydrogenated vaeaney, whieh is also predieted to be perfeetly 
three-fold symmetrie .1221 However, in reeent STM experiments and DFT simulations of Ar+ ereated 
hydrogenated vaeaneies on graphite, several different eonfigurations of hydrogen are shown with 
varying image eontrast,l2Hl where we see remarkable eonsisteney in the appearanee of the bright, 
triangular defeets. These eonsiderations suggest that the bright defeets we observe are indeed 
single vaeaneies. 

Briefly returning to the dark defeets, the white arrows in Fig. 3f indieate that these defeets do 
not disturb the graphene lattiee, whieh is eontinuous over the defeets. We also notiee that many 
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of these types of defeets appear to be swept to the edges of graphene islands during growth (see 
Ref.EUand Fig. la, b, and e). These observations, along with the similarity between these defeets 
and those on the bare Cu areas, and the seattering of the Cu surfaee state eleetrons (Fig. 2e, 3f) 
are strong evidenee that these are defeets in the underlying Cu lattiee.l^ The larger dark area in 
the eenter of Fig. 3f appears to be quite different, and may be a defeet eluster in the graphene. 
The upper inset of the figure shows this region with a Laplaeian filter applied to emphasize the 
surrounding graphene lattiee, whieh is elearly deformed at the edges of this feature. 

Lastly, we performed tunneling speetroseopy to probe how the defeets influenee the eleetronie 
strueture of the graphene (Fig. 4). Point speetra taken on the eenters of the dark and bright defeets 
in Fig. 4a are eompared to pristine graphene in Fig. 4b. The primary feature in these speetra is the 
band edge of the Shoekley surfaee state in the underlying Cu(l 11), whieh is shifted to -0.31 V for 
graphene-eovered Cu.EU There are slight differenees in the sharpness and magnitude of the band 
edge in dl/dV between the speetra on and off the defeets. These differenees are emphasized in a 
spatial map of the dl/dV signal at -0.31 V (Fig. 4e). The bright defeets show a greatly enhaneed 
LDOS that extends symmetrieally around the defeet with a FWHM of ~3 nm. The dark defeets 
have a smaller enhaneement of LDOS at their eenter, then a diminished LDOS in a ring shape 
surrounding the defeet with a diameter of --1.2 nm. Ring-like struetures in dl/dV maps like these 
have previously been attributed to ionization rings of eharged defeets, sueh as Co on graphene!®! 
and Mn dopants in GaAs. I®! Additional studies, sueh as a bias dependenee of the ring diameter, 
are needed to test whether this ring strueture is due to ionization, whieh would be surprising given 
that the defeets appear to be strongly eoupled to the metal substrate. Figures 4d-e show eolor 
maps of dl/dV speetra, taken point-by-point by moving the tip along the lines indieated in Fig. 
4a. These data indieate that the enhaneed band edge LDOS extends only ~lnm from both types 
of defeets.There are no other prominent defeet-speeifie eleetronie states within the range probed 
here. 

Tunneling speetroseopy of defeets in small, regularly shaped islands is further eomplieated 
by quantum eonfinement of the Shoekley 2DEG. Figure 5 shows one sueh island, with a step in 
the underlying Cu running through the eenter. This island eontains three bright defeets, elearly 
resolved in Fig. 5b under imaging at IV, but almost indiseernible in STM imaging at 0.05 V (Fig. 
5e). Tunneling speetra taken in a line over the upper bright defeet in Fig. 5b again shows the 
Shoekley band edge at -0.31 V as the primary feature eommon to all the speetra (Fig. 5d). In 
these data, however, sharp peaks in the LDOS exhibiting strong spatial and voltage dependenee 
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are also prominent (red features in Fig. 5d). These peaks arise due to quantum eonfinement of 
the surfaee state eleetrons within the small and regularly shaped island, similar to those of the 
quantum eorral.El Again, a elear signature in the speetroseopy due to the bright defeet is notably 
absent. For referenee, the line profile showing the bright defeet apparent height at Vbias =1 V is 
overlaid on the speetra in Fig. 5d. Any modifieation to the LDOS due to the defeet is diffieult to 
distinguish from that of the quantum-eonfined Shoekley eleetrons. 


IV. CONCLUSION 

We have reported a detailed STM study of natively oeeurring defeets in UHV grown graphene 
on Cu(lll). We foeused on two types of point defeets we eonsistently observe within otherwise 
pristine graphene islands and elassify them as bright and dark, aeeording to how they image at 
most biases, but speeifieally at V^ias =1 V. Atomie resolution images aehieved with funetional- 
ized tips reveal that bright defeets have a trigonal symmetry surrounding a C site in the graphene 
lattiee, whieh is rotated by 60° depending on the sublattiee assignment of the C site. These de¬ 
feets are identified as single vaeaneies arising in the as-grown graphene. The graphene lattiee is 
undisturbed by dark defeets, suggesting that these are point defeets in the underlying Cu lattiee. 
dl/dV speetroseopy of the defeets shows their primary eleetronie signature is in modifying the 
Shoekley surfaee state of the Cu(lll) substrate. More subtle effeets of the bright lattiee defeets 
on graphene’s LDOS may not be observable for this experimental system. On the other hand, 
experiments on UHV grown graphene have the advantage of eliminating possible environmental 
effeets, where eontaminants eould easily funetionalize C vaeaneies. This experiment also provides 
a study of graphene point defeets arising from a growth proeedure rather than ion bombardment 
or irradiation. We expeet this work to provide new perspeetive on the ereation and imaging of 
point defeets in graphene and motivation for pursuing growth as a means to tune defeet density for 
eontrollable doping and funetionalization for graphene deviees. 
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Figure 1: STM images of UHV-grown graphene islands on Cu(lll). Yellow arrows indicate examples of 
bright defects and blue arrows indicate examples of dark defects. The red arrow in (b) points to a cluster of 
dark defects and the green arrow in (e) points to a graphene grain boundary. Images in (b-e) are topographic 
images overlaid with images treated with a Laplacian filter to bring out small differences in height. All scale 
bars are 10 nm. Images (a-d) taken at 14ias =1 V and I =0.2 nA. Image (e) taken at 14ias =1 V and 1=1 nA. 


13 

















Figure 2: Graphene point defects, a) STM image at 1 V, 0.2 nA of a graphene island with 11 bright defects 
and 3 dark defects b) and c) depict example apparent height and apparent depth profiles, respectively, of the 
defects marked with horizontal white lines in (a) The vertical grey bar in (b) represents the standard devi¬ 
ation around the mean apparent height for bright defects in (a), d-f) Bias dependence of defect appearance 
for images taken at 0.2 nA and d) IV, e) 50 mV, and f) -IV for the region outlined by the yellow box in (a). 
Color scale varies from 0 (black) to 36 pm (yellow). 
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Figure 3: Atomic resolution imaging of point defects, a-c) STM images of a two bright defects (blue dashed 
circles) imaged at 14ias =10 mV and I =\ nA with a) a metal-terminated tip (Tip A), b) an unknown non- 
metal terminated tip (Tip B), and c) a CO terminated tip (Tip C). (Scale bar is 2 nm) d-e) Same defects as 
(a-c) imaged with Tip B using tunneling conditions that produce very high resolution of the defect structure 
and surrounding lattice. The lattice overlay shows that the defect is centered on a C site. (Vbias =10 mV 
and I =0.2 nA, scale bar is 1 nm) f) Atomic resolution image with a non-metal terminated tip (Tip D) and 
lattice overlay of a large graphene region with a number of dark defects and three bright defects. Lower 
insets: Magnified view of the two orientations of triangular bright defects showing their corresponding 
sublattices (red and blue). Upper inset: Laplace filtered image showing lattice deformations in the central 
dark defective region. (Ubias =106 mV and I =0.4 nA, scale bar is 2 nm) 
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Figure 4: Spectroscopy of point defects, a) Topographic image of a graphene island with 5 bright defects 
and 3 dark defects (one on upper terrace). A 1.5 nm Moire pattern is visible. (Vbias =-0.31 V, I =0.5 nA) 
b) dl/dV point spectroscopy to measure the local density of states on dark (red) and bright (blue) defects 
compared to bare graphene (black), at locations marked in (a). Dashed black line at -0.31 V marks the 
Shockley surface state band edge, c) Spatial dependence of dl/dV of region in (a) at 14ias =-0.31 V and 
I =0.5 nA. d-e) Spatial and voltage dependence of dFdV in dense lines (representing 30 spectra each) over 
d) a dark defect and e) a bright defect in (a) (yellow dashed arrows). 
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Figure 5: Quantum confinement effects on spectroscopy, a) STM image of a small graphene island with 
regular hexagonal edges. A Cu step edge with a 120° bend bisects the island, b-c) STM image showing 
a magnified view of a section of the upper terrace of the island (yellow dashed rectangle in (a)) with three 
bright defects imaged at / =1 nA and b) Vbias =1 V and c) 50 mV. Standing waves due to quantum con¬ 
finement of Shockley surface state electrons is visible in (c). d) dl/dV spectroscopy in a line over a bright 
defect with topographic line profile at 1 V overlaid (pink dashed arrow in (b)). Quantum confinement peaks 
show up as spatially varying, bright peaks in dl/dV. Red arrow marks the Shockley surface state band edge, 
blue and green arrows mark the voltages at which images in (b) and (c) were taken, respectively. 
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